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Abstract 

Properties of the proton and neutron are studied in partially-quenched chiral 
perturbation theory at finite lattice spacing. Masses, magnetic moments, the 
matrix elements of isovector twist- 2 operators and axial-vector currents are 
examined at the one-loop level in a double expansion in the light-quark masses 
and the lattice spacing. This work will be useful in extrapolating the results of 
simulations using Wilson valence and sea quarks, as well as simulations using 
Wilson sea quarks and Ginsparg- Wilson valence quarks, to the continuum. 



June 2003 



1 



I. INTRODUCTION 



Impressive progress is currently being made in understanding properties and interactions 
of the low- lying hadrons using lattice QCD. However, computational limitations necessitate 
the use of quark masses, m q , that are significantly larger than those of nature, lattice spac- 
ings, a, that are not significantly smaller than the physical scale of interest, and lattice sizes, 
L, that are not significantly larger than the physical scale of interest. In order to make a 
connection between lattice QCD calculations of the foreseeable future and nature, extrap- 
olations in the quark masses, lattice spacing and lattice volume are required. Assuming a 
hierarchy of mass scales, 

L' 1 < m q < A x < a' 1 , (1) 

where A x is the scale of chiral symmetry breaking (a typical QCD scale), and working in 
the infinite volume limit, the appropriate tool for incorporating the light quark masses and 
the finite lattice spacing into hadronic observables is effective field theory (EFT). 

Chiral perturbation theory (x?T) provides a systematic description of low-energy QCD 
near the chiral limit and is therefore an EFT which exploits the hierarchy m q <C A x . 
This technology has been extended to describe both quenched QCD (QQCD) [1-5] with 
quenched chiral perturbation theory (Q;\PT) and partially-quenched QCD (PQQCD) [6-10] 
with partially-quenched chiral perturbation theory (PQxPT). It is hoped that future lattice 
simulations can be performed with sufficiently small quark masses where the chiral expan- 
sion is convergent, and can be used to extrapolate down to the quark masses of nature. 
Recently, meson and baryon properties have been studied extensively in both Qx?T [4,5] 
and PQxPT [11-14]. The effective field theory (EFT) describing the low-energy dynamics 
of two-nucleon systems in PQQCD has also been explored [15,16]. 

In order to construct the low-energy EFT at finite lattice spacing, written in terms of 
the hadronic fields, one must first construct the underlying lattice theory, written in terms 
of the quark and gluon fields. The lattice theory and the continuum theory coincide in the 
a = limit, but away from this limit the theories differ. The lattice breaks the Lorentz 
group (0(4) in Euclidean space) down to the discrete symmetry group of the lattice, which 
we will take to be the symmetry group H(4) of a hypercubic lattice. As first discussed by 
Symanzik [17], the strong- interaction Lagrange density at a ^ will receive contributions 
from an infinite series of operators, ~ J2 a k O^ +k \ Therefore, the contribution from terms 
of 0(a n ) to a given strong-interaction observable will, according to eq. (1), be suppressed by 
factors of ~ a™ A™. For Wilson fermions [18], where chiral symmetry is not a good symmetry, 
it is straightforward to show that at 0(a) the Symanzik Lagrange density has the form, once 
appropriate redefinitions and renormalizations have been performed, 

C = ^(p+m q )^j + ac s Jpa^G^ + ... (2) 

where c sw is the Sheikholeslami-Wohlert [19] coefficient that must be determined numerically. 
For lattice fermions that satisfy the Ginsparg- Wilson (GW) condition [20], such as Kaplan 
fermions [21] and overlap fermions [22], where chiral symmetry is a good symmetry, the 
coefficient of the Sheikholeslami-Wohlert [19] term vanishes, c sw = 0. The power-counting 
we will use in this work treats both m q and the lattice spacing a as small. The small 
dimensionless parameters that we will use in our low-energy EFT are 
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in 

p 2 ~ir ~ aA x ~ T2" > ( 3 ) 

where <9 represents the derivative operator. 

Following earlier work of Sharpe and Singleton [23] (see also Ref. [24]), Rupak and 
Shoresh [25] have extended \PT to 0(p A ) including the effects of a finite lattice spacing 
at 0(a) for Wilson fermions, and computed the Goldstone-boson masses. Together with 
Bar [26] these same authors have generalized the results to "mixed" actions where different 
types of lattice fermions are used for the sea and valence quarks. Recently, this work was 
extended to 0(a 2 ) for both Wilson and mixed actions [27,28]. The pion decay constant has 
also been computed to 0(a 2 ) for Wilson fermions in QCD by Aoki [28]. When considering 
the matrix elements of operators coupled to external sources, such as the axial-vector cur- 
rent matrix elements or the matrix elements of twist-2 operators, there are contributions 
at 0(a) from the operator itself, in addition to the Sheikholeslami-Wohlert term in the 
strong-interaction sector. In this work we consider the leading 0(a) corrections to nucleon 
properties. We compute the contributions to the nucleon masses at 0(p 3 ), to their magnetic 
moments at 0(p), to their isovector axial- vector matrix elements at 0(p 2 ) and to the matrix 
element of the n = 2 isovector twist-2 operator at 0(p 2 ). 

II. PQQCD AT FINITE LATTICE SPACING 

The Symanzik effective Lagrange density at 0(p) which describes the quark-sector of 
PQQCD for two light flavors is 

C = Q \ % Jp-m Q \ Q + Q o^Gp, A Q Q , (4) 

where the valence, sea, and ghost quarks are combined into the column vector 

Q = (u,d,j,l,u,ct) . (5) 



The u and d are valence quarks, the u and d are ghost quarks, and the j and / are sea quarks. 
The mass matrix, rriQ, is ttiq = di&g(m u ,md,mj,mi,m u ,md), and the Sheikholeslami- 
Wohlert coefficient matrix is Aq = a diag(c^, c^Q , c^, c^, c^Q, [25,26]. As men- 
tioned previously, when both the valence and sea quarks are Wilson fermions c^Q = c^J, 
but when the valence quarks are GW fermions while the sea quarks are Wilson fermions, 
C W = 

The graded equal-time commutation relations for two fields are 

Q?(x)Qf (y) - (-) Wfc Q£ f (y)Q?(x) = rt^ 3 (x - y) , (6) 

where a, (3 are spin-indices and i, k are flavor indices. The objects i]k correspond to the parity 
of the component of Qk, with r]k = +1 for k — 1, 2, 3, 4 and rjk = for k — 5,6, and the graded 
commutation relations for two Q's or two Q^s are analogous. The left- and right-handed 
quark fields, Ql,r in eq. (5), transform in the fundamental representation of SU(4:\2) L R , 
respectively. The ground floor of Ql transforms as a (4,1) of SU(4) qL <g> SU(2)q L while 
the first floor transforms as (1,2), and the right handed field Qr transforms analogously. 
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In the absence of the quark-mass and Sheikholeslami-Wohlert terms, itlq = Aq = 0, the 
Lagrange density in eq. (4) has a graded symmetry U(A\2)l <8> U(4\2)r, where the left- and 
right-handed quark fields transform as Ql — > UlQl and Qr — > UrQr, respectively. The 
strong anomaly reduces the symmetry of the theory to S'£/(4|2) L (g>,S , L r (4|2) R (g>£/(l)y [10]. It 
is assumed that this symmetry is spontaneously broken according to the pattern SU(4\2)l<S> 
SU(A\2) R <g> U(l) v -> S77(4|2)y <g) C/(l)y so that an identification with QCD can be made. 



A. The Pseudo-Goldstone Bosons 



In order to construct the Lagrange density describing the dynamics of the pseudo- 
Goldstone bosons at 0(p 2 ), we allow mg and Aq to transform under the graded chiral 
group [23,25,26]. This leads to 



tl r-2 j-2 

C = J —str f d^E^E 1 + A M ^str f m Q ^ + m Q E 1 + A A ^str 



AqE* + AqY, 



(7) 



where a$ and m are quantities that do not vanish in the chiral limit. The meson field is 
incorporated in E via 



cxp 



2JJ> 



M 

X 



M 



(8) 



where M and M are matrices containing bosonic mesons while x an d X^ are matrices con- 
taining fermionic mesons, with 

J° 

' M = ( 7 1 V L n J ) , x = ( Xvu A7rT AJ " AjLT ) • (9) 



/ Vu 



M 



IT 

7° 



7T 

J+ 
1° 



^7 



Vi I 



Vu 
7T~ 



7T 

Vd 



XtT + 



Xj- 



Xl+ 
Xl° 



The convention we use corresponds to / ~ 132 MeV, and the charge assignments have 
been made using an electromagnetic charge matrix of Q( p< ?) = |diag (2, — 1, 2, — 1, 2, — 1). 
The singlet field is defined to be $ — str ( $ ) / \/2, and its mass m can be taken to be of 
order the scale of chiral symmetry breaking, m — > A x [10]. Hence the parameters a$ and 
mo decouple from the low-energy EFT in PQxPT [10]. 

It is straightforward to show that the meson masses resulting from this Lagrange density 

are 



m 



ud 



m 



m 



uu 
2 



Am (m u + m d ) + 2\ A a c 
2 Am rn u + 2 \a a c 



(V) 
sw 



(V) 



JU 



m 



A M (rrij + mu) + \ A (a c\J + a c\J) 
Am (rrij +mi) + 2X A a c 



(S) 



(10) 



and so forth, where m 2 ab denotes the mass-squared of a meson containing a (anti-) quark 
of flavor a and one of flavor b (either valence, sea or ghost). The meson masses have been 
computed out to 0(m 2 a) in Refs. [25,26]. 
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B. The Nucleons and A-Resonances 



The free Lagrange density for the 70-dimensional baryon supermultiplet, B^ contain- 
ing the nucleon and for the 44-dimensional baryon supermultiplet, %^ k1 containing the 
A-resonances is [11,12], at leading order (0(p)), 

£ = i(Bv VB) + 2a£ Q) (BBM+) + 2[3 { m Q) (BM+B) + 2<r£ Q) (SB) str (M+) 
+ 2«i PQ) (BBA+) + 2(3 { A PQ) (BA+B) + 2a A PQ) (BB) str (A+) 

- i (T% ■ VT») + A (T%) - 2 7 £ Q) (T»M + T,) - 2<f£ Q) str (M + ) 

- 2 7 i PQ > (T»A + %) - 2*f» {T%) str (A) , (11) 

where A is the mass splitting between the 70 and 44, M.+ = \ {^ttiq^ + £ — V^, 

and A + = \ (p Aq^ + £Aq£^J . The terms that arise at 0(p) from the Sheikholeslami- 

Wohlert [19] operator have coefficients a ( A Q \ (3 ( A Q) , a ( A Q \ 1a pq \ and a ( A Q) . 

The Lagrange density describing the interactions of the 70 and 44 with the pseudo- 
Goldstone bosons at LO in the chiral expansion is [4], 

£ = 2a (BS^BA„) + 2/? (BSVAvB) + 2H (j'S^A^) 

+ sjhz [ (T U A U B) + (EA V T U ) ] , (12) 

where S* M is the covariant spin- vector [29-31]. Restricting ourselves to the valence sector, 
we can compare eq. (12) with the LO interaction Lagrange density of QCD, 

£ = 2g A NS^A^N + 9l NS^N tr [ A^ ] + g AN [ T bc ' u A d av N b t cd + h.c. " 

+ 2saa T V S»A^T U + 2g x T S" T v tr [ A^ ] , (13) 

and find that at tree-level 

4 1 2 1 

« = gfiU +3^1 » P = 3^1 - ^ 9a ' n = 9AA ' C = ~ 9AN ' ^ 

with = 0. Considering only the nucleons, and decomposing the Lagrange density in 
eq. (13) into the mass eigenstates of the isospin-symmetric limit, 7r ± , 7r° and i], we have 

_ - \/2 — 

C = 2g A NS»A,N + ^ (g A + 9l ) NS»N , (15) 

where A^ is the axial- vector field of pions only (excluding the isosinglet meson). In the 
isospin-symmetric limit, with the mass of the 77 being of order ~ A x , all expressions must 
be independent of the coupling g±. At the order we work to in this paper, higher-order 
interactions do not contribute. 
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III. NUCLEON MASSES 



The mass of the i-th baryon in the 70-dimensional baryon supermultiplet has an expan- 
sion in m q and a, of the form 

Mi = MM - iWfV) - Mf /2) (/i) + ... , (16) 

and we will be interested only in the proton and neutron masses, i.e. i — p,n. The super- 
script denotes the order in the expansion, i.e. M^ 2 \ji) denotes a contribution of 0(p 3 ). 
The term Mo(/i) is the same for all baryons in the supermultiplet, and is non-zero in the 
m,Q, a — > limits. The a = values of M- (//) and M- (/x) can be found in Ref. [12] for 
arbitrary quark masses. For a ^ we find that V) 

becomes 

M« = \m, ( 5a ™ + 2/?™) + \m d ( a ™ + 4/?™) + 2<r™ K + ™<) 

M« = im„ («™ + 4/3™) + im, (So™ + 2,^ Q >) + 2«x<f» (m, + m,) 

+ 2ac«? («<«»+ + . (17) 



At this order, the finite a contributions are the same for the proton and neutron. This has 
to be the case as the lattice spacing, being the same for the u and ci-quarks, transforms 
as an isoscalar. The M^ 2 \n) contributions, arising from one-loop diagrams in PQxPT, 
have exactly the same form as in eqs.(38) and (40) of Ref. [12], however, it is understood 
that the meson masses are evaluated with the relations in eq. (10), and hence have implicit 
dependence on the lattice spacing. We see that this introduces non-analytic dependence on 
the lattice spacing in the chiral limit. 

As the explicit expressions for the loop contributions are quite long, we will simply quote 
the results in the isospin limit (i.e. the sea quarks are degenerate, but different in mass from 
the degenerate valence quarks and ghosts), 
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M (3/2) = 1 9A 

P " 8tt/ 2 I 12 



,2 



9m 2 ss mvv + 16m|y — 7my V + ^- 9m z ss mvv + Wmg V — 19my V 

A. 



9a9i 
6 



29 



9m 2 ss m vv + 4m| y - 13m d vv + ^AK [F vv + F sv ] ) , (18) 



3tt 



where S denotes a sea quark and V denotes a valence quark. The function F c = F(m c , A, ji) 



is 



F(m, A, fi) = (m 2 - A 2 ) j VA 2 - m 2 log 





(19) 

where fj, is the renormalization scale, and the meson masses depend upon light-quark masses 
and the lattice spacing through eq. (10). Note that we use dimensional regularization with 
MS to regulate divergent integrals. Of course, all regulators must give equivalent results. 
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In the QCD and isospin limit, mj,mi,m u ,md — > m, while at finite a, the nucleon mass 



is 



m n = M - 2 m + P ( m Q) + 2-r } ) - 

-2ac2? + " 4 a eg,™ . 



3 2 3 4(7 AA r 



(20) 



IV. NUCLEON MAGNETIC MOMENTS 

The most general analysis of the nucleon magnetic moments requires us to determine the 
vector- current operator out to the order in the lattice spacing that we are working. While 
the O(a) corrections are known [32], in order to compute the magnetic moments at 0(p 1 ^ 2 ) 
only the continuum limit of the vector-current operator is required. Up to 0(p 1 ^ 2 ), it is 
convenient to write the magnetic moment of the i-th nucleon as 



Hi = OLi + 



M, 



N 



Pi + ft ] + 



(21) 



where the constants, a i: $ and f3[ in PQQCD are given in eqs.(49)-(51) in Ref. [12]. The 
only modification required at finite lattice spacing is to use eq. (10) in evaluating the meson 
masses. Therefore, we see that the lattice spacing first appears through the meson masses at 
order 0(p 1 ^ 2 ), and hence the leading contribution of the finite lattice spacing is non-analytic 
in the lattice spacing in the chiral limit. 
In the isospin limit one finds 



g 2 

~ [8m sv + m vv ] 



[m sv ~ m vv ] 



[m sv ~ m vv ] 

lo 



+ {0j + qi) [m S v ~ mvv] 



9A9\ 



9l 



9an 
9 



-1T VV + -(qj + qi) {Tvv - ?sv) 



Y [4msv + 5m vv ] + 9 ^ 91 [m sv - m vv ] - 

If 5 
+ (Qj + qi) [m S v - m vv ] - \2g 2 A + g A g 1 + -g\ 



2 9 f 



[msv ~ m vv ] 



9an 
9 



Tvv + ?sv + -(qj + qi) {Fvv - Fsv) 



(22) 



where the function Ti = T{m^ A, ji) is 



7rjF(m, A, //) = V A 2 — m 2 log 




(23) 



We have used Q = diag (+§, — |, qj, qi, qj, qij, for the electromagnetic charge matrix in 
PQQCD [11,12]. The charges of the sea quarks and ghosts, qj and qi, arise due to the 
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fact that their charge assignments are not unique, only constrained by the requirement that 
electromagnetic observables computed in PQQCD reproduce those of QCD in the QCD 
limit [11,12,36]. It is clear from the expressions in eq. (22) that this is indeed the case. 

In the isospin-symmetric QCD limit at finite lattice spacing, the magnetic moments 
become [33,34] 



\h = A*o + Hi 



= fJ>o - A*l + 



M 



N 



4tt/ 2 
M N 



2 _i_ 2 f 

9 A m K ' g 9 AN •> ?r 

2 ^2 — 

9A m 7T + g fi'AAf -F K 



(24) 



The isoscalar and isovector magnetic moment contributions from the LO dimension-5 oper- 
ators are /i and fj,i, respectively, and are independent of the lattice spacing. 



V. NUCLEON AXIAL MATRIX ELEMENTS 

The leading effects of a finite lattice spacing on the matrix elements of the axial-vector 
current enter at 0(p 2 ) from both one-loop diagrams and from local counterterms. In addition 
to the contribution from the Sheikholeslami-Wohlert term, there is also a contribution from 
the 0(a) corrections to the axial-current operator. As discussed in detail in Ref. [35], there 
are only two operator structures that contribute to the axial-current operator at 0(a), using 
the notation of Ref. [35] , 

Oj,, = 1 r a l5 (t 5„) q , = q r a W m q q , (25) 

in QCD, where D^= D ^ — D^. The appearance of m q in 0%^ renders it 0{p A ) in the 
power-counting, and so we neglect it, leaving the contribution from Oj^. The extension to 
PQQCD is straightforward. In analogy with the electromagnetic interaction, the isovector 
axial-charge matrix must be extended to PQQCD and the axial charges of the sea quarks 
and ghost must be defined as the extension is not unique [36,11,12]. Requiring the axial 
charge matrix to be supertraceless implies that the most general extension to the leading 
operator is, e.g., 

t 3 ^r 3 = diag (1,-1, yj,yi, yj,yi) . (26) 

In order to construct the nucleon axial matrix element at leading order one uses the spurion 
construction in which T a L — > Lr a L L^ and T a R — > Rr a L R) , where the axial-current operator is 
decomposed into contributions from the left- and right-handed quark fields. 

At subleading order, 0(p 2 ), the contribution to the axial current in PQQCD is 

6AI = <3t% 7 75 (iD^Q , (27) 

where, for example, the matrix t\ 7 is 

t\ 7 = a diag (c% ] , -c% ] , yjcffl , yic^j , y^ , yic% } ) , (28) 
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for the charge matrix in eq. (26), where c^-J and the coefficients of the axial-current 

corrections for the sea quarks and valence quarks, respectively. Under chiral transformations 
the spurion field transforms as t\ 7 — > Lr\ 7 W, like the quark mass matrix. 

The leading-order contribution to the matrix elements of the axial current arises from 
the operators [11,12] 



( p °)j; 5 = 2a (BS,B r« + ) + 2(3 (BS, + 2H (T V S» t%+T v ) 



+ 



(29) 



where t£ + = \ ( £T a ^ + ^r a £ ) . 

At 0(p 2 ) there are three different contributions, that we denote as ™8j$ a , {PQ) Sj { ^ a 
and ( PQS> 5j^l,' a . The contribution from a single insertion of the leading 0(a) corrections to 
the axial-current operator is 

{PQ WS' a = 27ai (BS,B rX^ + ) + 2 7A2 (BS, t\^ + B 



(30) 

where t A j^ + — \( ^aj^, + &Ta,7& ) • The remaining two contributions arise from a single 
insertion of the light-quark mass matrix and from the Sheikholeslami-Wohlert term, which 
are of the form 

(PQ) & 3) ' a = 2 [ &i,r B k]l { t| + , T + S,B njk 

+ 6 2 ,r (-)^+*)^+^) { rl + , r + }£ S„B ijn 

+ v B fc - ((r« + ); (r +) ; + (r + )J 

+ 6 5 ,r R^+^ B~ kji (rf + Y. (T + )? + 6 6 , r fe + V S^ fe str (r + ) 

+ b 7 , r (r» + )^ str (r+) 



+ 6 8 , r IT* S^-* str (rf + T + ) 



(31) 



where T = Ai and ^4. 

We write the axial matrix elements as [12] 



(N b \™j,, 5 \N a ) 



Pah + 



167T 2 / 2 



Vab 



Pab^Z [ W a + W b } + Vj T] ( £ + Vl 77, 



ah 



+ c ab + y c% + Vl c2 + d ab + y 3 <%> + yi d a l > 2 U b S,U a , (32) 



00 



where the constants p ab , r] ab , r]^ h \ w a , c ab , c a K b \ can be found in Ref. [12] with bj — > bj yM . In 
extending the r + isovector operator from QCD to PQQCD, one can simply replace the r 3 in 
the upper 2x2 block of r 3 with r + , as described in Ref. [12]. The meson masses in p ab) rj ab , 
7] a k \ w a are understood to be evaluated at finite lattice spacing via eq. (10). The constants 
d ab and d$ arise from ^8j^ a and and are 



,(*) 



(V) 

+ ( 2 7A)1 - 7Ai2 ) , 

= d (i) = d U) = d (i) = b ( (S) _ (V)\ 

u pp u pp u nn u nn U S,A \ u-^ S w uu sw J ) 

dnp d nn dpp , <^np ^np ^ " (<^3) 

In the isospin-symmetric QCD limit at finite lattice spacing, the proton isovector axial 
matrix element is 

(pl&b) = 9 a - (^9 a (l + 2^) + (2# A + ^aa)^jv -4 - ^9A9\ N K^j 

+ y (-26 1;M + 46 2jM - 6 3jM + 6 4jM + 26 5jM - 26 6jM + 46 7jM ) 

ac^) 2ac (s) 
+_^L [ _ 2 & M + 4b 2A - h A + 6 4 ,a + 2b 5A ] + — [2b 7A - 6 6 , A ] 

+ °^f- ( 2 7A;1 - 7A, 2 ) + (%• + yi)bs,A (acg? - acg?) , (34) 
where L w = m 2 log (m 2 //i 2 ), J,,- = ./(m,,-, A, /i) is 

J(m, A, „) = (m 2 - 2A 2 ) log (^) + 2Av^^ log ~ " ^ + , (35) 

v ' \fi z J \A + V A 2 -m 2 + tej 

and K n = K(m n , A, //) is 

2 a 2 A, /™ 2 \ 2 . - /A - VA 2 -m 2 + ie\ 



K(m,A,/i) = (m 2 --A 2 )log(^j + -AVA 2 -m 2 log 



,/i 2 / 3 VA + v / A 2 -m 2 + k/ 

2m V /Al 2i A - VA 2 - m 2 + 

+ 3Al mWA l0g ( A + VA 2 -^ + W J ' (36) 

where fi is the renormalization scale. The last contribution in eq. (34) is somewhat interesting 
as it depends upon the extension from QCD to PQQCD and also upon the choice of lattice 
fermions and does not vanish in the isospin-symmetric QCD limit. 



VI. MATRIX ELEMENTS OF THE ISOVECTOR TWIST-2 OPERATORS 

The forward matrix elements of twist-2 operators play an important role in hadronic 
structure as they are directly related to the moments of the parton distribution functions. 
In QCD the long-distance contributions to these matrix elements have been computed order- 
by-order in the chiral expansion using %PT [37-40] and have been applied to results from 
both quenched and unquenched lattice data [37], with interesting results. Further, the 
analogous contributions in QxPT and PQ%PT have been computed in Refs. [11,12,41]. 

In QCD, in the continuum limit the twist-2 operators are 
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O^L ..., n = 9^ 7^ (i B M2 ) ... (i B M „ ) g - traces , (37) 



where it is understood that the operator is symmetrized on its Lorentz indices. For the for- 
ward matrix elements in the proton and neutron, as are relevant for deep inelastic scattering, 
we only need to consider 6 = 3. The extension to PQQCD is straightforward, 

{PQ) { ^\ = QT b 7/ll (i B, 2 ) ... [i D, n ) Q - traces , (38) 

where r b has the same form as that in eq. (26), and we will use the same charges, but it 
should be remembered that they are unrelated to those of the axial-current operators. 

The higher-dimensional operators that enter at finite lattice spacing are, in general, quite 
complicated. This is due to that fact that the operators must be classified not only by charge 
conjugation, parity and so forth, but also by the representation theory of the hypercubic 
group, H(A). Such a classification has been performed up to n = 4 [42], but as n increases 
the number and complexity of the operator basis increases significantly. 

The n = 1 operator, Off' 3 , is the isovector vector-current operator with 0(a) correc- 
tions [32] 

q r 3 ^ m q q , d v ( q r 3 o» v q) . (39) 

The forward matrix element of the second operator vanishes for obvious reasons and the 
first operator is 0(p A ) in the expansion, and thus there are no operator corrections to Off' 3 

to the order we are working. Therefore, there are no modifications to the matrix element of 

n(i)>3 

By contrast, the forward matrix element of the n = 2 operator, 

Off' 3 = qr 3 ^ (i 5„) q , (40) 
does receive corrections at 0(p 2 ) [32]. The correction to Off' 3 at 0(p 2 ) in QCD is 



SO$* = acff q r 3 a 



i Dv , i Dx 



q + acff qr 3 {iD^ , iDv} q ■ (41) 
When extended to PQQCD, this correction becomes 

{PQ) SOff' 3 = Q r 3 A1 <j, x [ i D v , i Dx ] Q + Q r% 2 {tD,,iD v }Q , (42) 
where the charge matrices are 

Tj,, = « dtag(<f -<f ><">, y, cf<»„ cp>< s u cf c?><"> ) , 

Tiu-.dhg^.-^.w 4 2,(S, , ! /, 4 2)(s) ,» £ m ,v, 4 2,, "» ) , (43) 

for the charge matrix given in eq. (26). The coefficients cf^^ 1 and c^ S ^ are the valence- 
and sea-quark coefficients, respectively. In order to construct the nucleon matrix elements 
we define the fields 
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T 



(2),3 



\ ( + ) , r?g + = \ ( ^ >2 £ + ?t\ 2 Z 



(44) 



In QCD, the matrix elements of the continuum operator Off' 3 are reproduced at leading 
order by 



„ . VflVu Nt* + N + 7 w ^ T" r| + T a 
— traces , 



,( 2 ) 



2 



T u r| + T M 



(45) 



and in PQQCD this becomes 



PQ 0^ - «?> 

+ 4 2) J 



B B r| + ) + AS 2) v>* ( B rf + B ) + 

( ^ i+ ^) + ( ^ i + ^) 



(2) 

7o ^ 



— traces 



(46) 



The 0(a) corrections to the operator, give rise to 0(p 2 ) corrections to the nucleon matrix 
elements of the form 



PQ SO W)* = «g> v,„ ( B B r?g + ) + /?g> v»v„ ( B B 

+ ag> v,v v ( B B r?g + ) + /?g> v,v u ( B rgg + B ) - traces , (47) 

while the contributions from a single insertion of the light-quark mass matrix and from the 
Sheikholeslami-Wohlert term at 0(p 2 ) are 



PQ 5(9 r,(2),3 = 2 



&g> a** { r« + , r + B njk 



+ 4 2 r R^'^ B*' (r^ + )J (r + )? B lnk 



-t- o 4 r 



^ (r +) ; + ( r + )J a, 

+ (-)*(*+*> £^ (r« + )' (r + )? £U + 6$ (r^ B ljk str (r + ) 



(ui+UiX'Jfc-HM) i3 fcji (r? + ) n %„ str(r+) 



+ b {2) 
-+- o 7 T 



+ 6$ B hjt B ijk str (rf + r+) 



f M f — traces , 



(48) 



where r = M and A 

The forward matrix elements between nucleon states can be written as [12] 



,3\. 



T 2 f 2 V " 



,(2)J 



167T 2 / 



(2),/ \ 

- p- 'Wi + yj rjl - y t r)l ) 



Ui v^v u Ui — traces , 



(2),i 



(49) 



where expressions for pf \ r/^' , Wi, r}^' k , cf^'° and the cf^ k can be found in Ref. [12] with 
the understanding that the meson masses are evaluated at finite lattice spacing according 
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to eq. (10). The constants df ] '° and df' k arise from the operators PQ 50°>^> 3 , PQ 50^' (2 ^ 3 
and PQ 50fy^' 3 , and are found to be 



d {2),0 = aC l 



(2)(V) 



(2) 
A 



(50) 



and those for the neutron are related by dffl' = —d^'° and = dfflj = dffl' 1 . 

In the isospin-symmetric QCD limit at finite lattice spacing these expressions reduce to 



( PQ o$>\ 



M ± V 

8tt 2 / 2 



Lit I - 



9an 
4tt 2 / 2 



5 (2) _ A^(2) 



27 



-CP 



+ o V + 40 2,M °3,M + °A,M + Z0 5,M Z0 6,M + 40 7,M J 



(2)(V) 



ac^ V) 



( 2 J 2) - 3 {2) \ 
+ \„JY) ( -2h {2) + Ah {2) - h {2) + ft (2 ' "™ " 2 

+ g aC sw ^ Z0 1,A + 40 2,A °3,A + °4, 

+ (%• + j/,) ( ac£J - ac£? )4 2 i 



r(2) , 2fi (2) 



3 

traces 



ac 



(S) 



(2) 
6,A 



26S 



(51) 



Like the matrix elements of the axial current, the choice of charges in the sea quark and 
ghost sectors remain at finite lattice spacing through the last term in eq. (51) with coefficient 

l(2) 
°8,A- 

VII. CONCLUSIONS 



As lattice QCD moves closer to its ultimate goal of computing strong-interaction ob- 
servables directly from first principles, effective field theory calculations must be developed 
in parallel in order to facilitate comparison with data and to make rigorous predictions. 
Significant attention has been paid to the chiral extrapolation of existing quenched and un- 
quenched data with the goal of making a connection between lattice calculations performed 
at unphysically large quark masses and nature. With this program maturing pleasantly, the 
time is ripe to address other extrapolations that need to be performed in order to make a 
rigorous connection with data: the continuum extrapolation, a — > 0, and the infinite volume 
extrapolation. 

In this work we have computed the leading effects of a finite lattice spacing, at 0(a), on 
some nucleon properties. One source of this dependence is the leading 0(a) corrections to the 
strong-interaction Lagrange density, i.e. the Sheikholeslami-Wohlert term. However, when 
considering matrix elements of operators, there are additional contributions from the 0(a) 
corrections to the operators themselves. If the lattice calculations are performed with lattice 
fermions that respect chiral symmetry then all of the finite lattice spacing contributions we 
have computed in this work will vanish. This will also be the case for C(a)-improved lattice 
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simulations. The continuum extrapolation of unimproved simulations of nucleon properties 
with Wilson quarks or mixed quarks, e.g. Ginsparg- Wilson valence quarks and Wilson sea 
quarks, can be performed with the expressions we have determined in this work. 
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